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Abstract. 

The time-dependent Dirac equation is solved using the three-dimensional Finite 
Difference-Time Domain (FDTD) method. The dynamics of the electron wave packet 
in a vector potential is studied in the arrangements associated with the Aharonov- 
Bohm effect. The solution of the Dirac equation showed a change in the velocity of 
the electron wave packet even in a region where no forces acted on the electron. The 
solution of the Dirac equation agreed with the prediction of classical dynamics under 
the assumption that the dynamics was defined by the conservation of generalized or 
canonical momentum. It was also shown that in the case when the magnetic field was 
not zero, the conservation of generalized or canonical momentum was equivalent to the 
action of the Lorentz force. 
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1. Introduction 

In our previous papers [TJ |2], the Finite Difference- Time Domain (FDTD) method, 
originally introduced by Kane Yee [3] to solve Maxwell's equations, was for the first 
time applied to solve the three-dimensional Dirac equation. The Zitterbewegung and 
the dynamics of a well-localized electron were used as examples of FDTD applied to 
the case of free electrons. The motion of electron wave packets inside and scattering 
from the potential step barrier or linearly dependent potential, arrangements associated 
with the Klein paradox j3], were used as examples of interaction with a scalar potential. 
In this paper, a FDTD study of the dynamics of an electron wave packet under the 
influence of a vector potential is presented. Such a dynamic behavior is most often 
associated with the Aharonov-Bohm effect [3]. 

As a manifestation of quantum mechanics, charged particles passing around a long 
solenoid can feel a magnetic flux even when all the fields in the region through which 
the particles travel are zero. The shifts in the phase of the wave functions describing 
the particles have been experimentally verified by its effect on the interference fringes 
P El El E] • Since, classically, there are no forces acting on the charged particles in the zero 
field region, the theoretical description of the Aharonov-Bohm effect contains a number 
of assumptions. They include assumptions on nonlocal features of quantum mechanics, 
the physical meaning of the vector potential, topological effects, etc. Generally accepted 
physical understanding is still lacking [lOj. While there is still an open question on the 
presence of classical forces responsible for the Aharonov-Bohm effect [TUl [EU H2] , on a 
macroscopic level they have not been observed |13j . 

Proper quantum-mechanical description of the dynamics of a relativistic charged 
particle involves the solution of the Dirac equation in the time domain. In addition 
to initial conditions, such a dynamics is defined only by the configuration of the 
scalar and vector potentials, and does not involve knowledge or any assumption on 
"classical forces". The solutions of the time-dependent Dirac equation, some of which 
are presented in this paper, can shed light and fill critical knowledge gaps on the 
theoretical and experimental interpretations of the mechanism of the Aharonov-Bohm 
effect, including the existence or non-existence of classical forces. 

2. Time-dependent solution of the Dirac equation 

The FDTD solutions of the time-dependent Dirac equation were obtained for the case 
when the electromagnetic field described by the four-potential A^ = {A (x), A(x)} was 
minimally coupled to the particle [HI [TS] 




(1) 



where 



H 



free 



—icha ■ V + j3mc 2 , 



(2) 
(3) 



int 



—ea ■ A + eA t 



FDTD solution of the Dirac equation and the Aharonov-Bohm effect 



3 



(4) 



and 
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The matrices a and (3 were expressed using 2x2 Pauli matrices a s and the 2x2 unit 
matrix I. 

In the FDTD method, the time dependent solution of the Dirac equation was 
obtained using updating difference equations [Tj. As an example, the values of ^1 at 
the position (iAx, jAy, kAz) and at the time step (n + 1/2) At were obtained using the 
equation 



1 2 (j J K)- 1 ( '" { > L - LK \u 



cAt 
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2AxC n (I, J, A) 
+ + 1, J, K) - - 1, J, K) - i(*J(J, J + 1, AT) 
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- i^(J, J, J, A) + J, A)^(J, J, A)], (5) 

where C n (J, J, A) = 1 + i§[mc 2 + eA%(I, J, A)]. The space and time were discretized 
using uniform rectangular lattices of size Ax, Ay and Az, and uniform time increment 
At. While it is not generally required, in the Eq. (|5| Ax = Ay = Az. Updating 
equations for ^2, ^3, and ^4 were constructed in a similar way. 

The dynamics of a Dirac electron can now be studied in any environment described 
by a four-potential N 1 regardless of its complexity and time dependency. In this paper 
we studied the dynamics defined by the vector potential A(x) 7^ 0. 

The Dirac equation is a differential equation of the first order and linear in d/dt. 
As in the case of Maxwell's equations, the entire dynamics of the electron is defined, 
only by its initial wave function. The dynamics of a wave packet used in this paper was 
defined by the initial wave function of the form 

/ 1 \ 



tf(£,0) = N 



E + mc 2 
2E 



V 





P3C 
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/ 



where N = [(2tt] 



3/2^.31-1/2 ; 



is a normalizing constant. Eq. m\ represents a wave packet 



whose initial probability distribution is of a normalized Gaussian shape. Its size is 
defined by the constant x , its spin is pointed along the z-axis, and its motion is defined 
by the values of the momenta Pi,P2, and p^. Some consequences of the initial localization 
of the wave packet on the overall dynamics of the electron were studied in Ref . [Tj [16] . 
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3. Validation of the computation: the dynamics of a wave packet in a 
strong uniform magnetic field 

The dynamics of a wave packet is very complex. The dynamics of a particle described 
by the wave packet in Eq. ^ depends on its localization, defined by the Gaussian 
component of the wave function, and its initial momentum, part of the wave function's 
phase. While an extensive study was done on the dynamics of the wave packet related 
to the scalar component A (x) of a four-potential [2], such a study does not validate the 
dynamics of the wave packet related to the vector component of a four-potential A(x). 

Applying the FDTD method to study the dynamics of a wave packet in a vector 
potential associated with a strong uniform magnetic field is not difficult. Classically, the 
dynamics consist of uniform rotational motion. In the relativistic quantum-mechanical 
description, however, even such a simple dynamics can validate the FDTD computation 
only up to a certain level. As pointed out in Ref. (17J "the dynamics is particularly 
rich and not adequately described by semiclassical approximations" . In Ref. [TS] it was 
demonstrated that in the presence of an external magnetic field the wave packet splits 
into two parts which rotate with different cyclotron frequencies, and after a few periods, 
the motion acquires irregular character. As a result of such a complex dynamics, when 
comparing the results obtained by the FDTD method and the computation described in 
Ref. p2] and [18], we could not expect more than a qualitative agreement. In addition 
to this qualitative agreement, of equal importance to the validation of the computation 
is the consistency of the results obtained for different vector potential gauges. 

In this paper, the motion of the wave packet described by Eq. ^ was studied in a 
uniform magnetic field oriented along the y-axis 



For this field the corresponding vector potential in a rotationally invariant gauge was 



The dynamics of the wave packet was obtained by solving the Dirac equation for this 
vector potential. 

While the probability densities \^\ 2 were calculated for the entire computational 
volume and at every time step, their values are shown here only on two planes, the 
horizontal or plane of classical particle motion, and the plane vertical to the plane 
of motion. The schematics of the planes relative to the wave packet motion and the 
orientation of the magnetic field are shown in Fig [TJ The position and the shape of the 
wave packet in the horizontal plane, as it moves along its first orbit, is shown in Figure 
[2j The initial position of the wave packet was at the center of the vector potential. Its 
initial momentum was pi = 0.53 MeV/c, making the motion relativistic. In order to 
force a relativistic electron to complete a full circle in the available computational space, 
the magnitude of the magnetic field was B$ = 10 s T. The field of such a strength is 
associated with the fields at the surface of the neutron stars. The classical orbit of the 



B = (0,50,0). 



(7) 



A = ^(-z,0,x). 



(8) 
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Figure 1. Schematics of the two planes used in displaying the probability densities. 
They are shown relative to the classical orbit of the electron and the orientation of the 
uniform magnetic field. 

electron in this field is r c i ass = p±/ {eBo) = 1.76 x 1CT 2 nm. As shown in Figure [2j during 
the first rotation the wave packet generally follows the classical orbit and disperses at 
the same time. The position of the center of probability of the wave packet relative to 
the classical orbit is shown in Fig [3} 

The complexity of the dynamics of the wave packet increased in the later stages of 
the motion. While the wave packet followed the circular motion, the probability density 
|^| 2 at some times assumed the spiral shape shown in Figure 4 increased or reduced its 
length, changed its rotational motion, and translated from one place to another. Overall, 
the characteristics of this dynamics were similar to the characteristics described in Ref. 
[T7] and |l8j . While the figures show some of the complex shapes of the probability 
density |^/| 2 , the richness of the electron motion can be better appreciated through the 
animation of the dynamics accessible on-line |19j . 

There was no intention of this validation to study the property of the confinement 
of the electron by the electrostatic potential. But, by using a scalar potential A Q at the 
boundaries of the computational volume to limit the electron motion in the direction of 
the magnetic field, due to the wave packet dispersion, we have inadvertently also solved 
the Dirac equation for the confined electron. We reported [2] that, contrary to the 
claim of the Klein paradox, there was no penetration of the wave packet into a potential 
barrier of a supercritical potential of height V satisfying the condition eV > E + mc 2 . 
The time evaluation of the solution of the Dirac equation shown in Figure [5] confirms 
this result. Figure [5] shows that the electron wave function, as the wave packet disperses 
and reaches a non-penetrable electrostatic potential, creates a bound state. Reduction 
of number of nodes in a bound state solution could be interpreted as a transition of the 
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particle from the higher energy state into the lover energy state as it orbits in a uniform 
magnetic field. 

Several additional tests of the FDTD method were performed using the dynamics 
of the electron in the uniform magnetic field. As expected, due to the normal spin and 
the magnetic field orientation, reversal of the orientation of the magnetic field resulted 
in the reversal of the direction of the rotation of the wave packet, keeping the same 
properties of the dynamics of motion. 

Of particular interest was testing the effects of the choice of gauge. The motion of 
the same wave packet was also studied for the uniform magnetic field oriented along the 
y-axis defined by the translationally invariant gauges 

A = B (—z, 0, 0), (9) 

or 

A = B o (0,0,x). (10) 

In both cases the wave packet persisted in a circular motion following the classical 
orbit. While the dynamics of the wave packet behaved as expected, using translationally 
invariant gauges has an additional importance on validating the FDTD computation. 
As seen in Eq. (J5l), the gauge in Eq. ^ couples to the \l/4 component and the gauge 



in Eq. (10) couples to the ^3 component. Similarly, the cross coupling exists for other 
components of the wave function 1 i r . As a result, the same dynamics should be obtained 
by different combinatorics of the components of the vector potential and the components 
of the wave function. This enabled for testing of possible inconsistencies in the FDTD 
updating equations. No inconsistencies were found. 

The dynamics of the wave packet motion in a uniform magnetic field was used here 
only as a validation of the FDTD method when a vector potential was applied in the 
Dirac equation. The same complexity of the quantum dynamics of motion as in previous 
publications [T71 [18] was shown. While used here only for computational validation of 
the FDTD method, this dynamics could be studied as a separate problem in more detail 
in the future. 

Finally, the complexity of the quantum dynamics of the wave packet in a uniform 
magnetic field studied here for three choices of gauge could be fully appreciated only by 
downloading the related animations [T9] . 



4. Wave packet dynamics in a vector potential created by two infinite 
solenoids 

The goal of this work is to study the dynamics of an electron wave packet under 
the influence of a vector potential associated with the Aharonov-Bohm effect [5]. 
Particularly, in this paper we present the dynamics of a wave packet obtained from the 
solution of the Dirac equation with a vector potential created by two infinite solenoids. 




Figure 2. The shapes and the positions of the wave packet in the horizontal plane 
during the first rotation. The circle represents the classical orbit. The animation can 
be accessed on-line [19] . 
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Figure 3. Positions of the center of probability of the wave packet in the horizontal 
plane during the first rotation. The circle represents the classical orbit. The animation 
of the dynamics can be accessed on-line [19] . 



The vector potential of a single infinite solenoid oriented along the y-axis can be 
written as 



A 



2ttR% 
2irr 2 



(-z,Q,x) 
'-z,0,x) 



for r < R 
for r > R 



where $ = BqtcRq is the magnetic flux inside the solenoid, Bq defines the strength of 
the magnetic field, R is the radius of the solenoid, and r = y/x 2 + z 2 is the distance 
from the center of the solenoid in the x-z plane. Outside two parallel infinite solenoids 
separated by a distance 2a and with opposite magnetic field orientation, the vector 
potential is then 

$ / z + a 
2^ \x 2 + (z + a) 



A 



x z + (z 



0. 



x 



X 



x z + [z 



X' 



(12) 



An example of the shape of this potential is shown in Figure [6j Outside the solenoids, 
the associated magnetic fields B = V x A of the vector potentials described by Eqs. 



(11) and (12) are zero. 



The time-dependent Dirac equation is now solved for the scalar component of the 
four-potential A = and the vector component A defined by Eq. p2| ). With this 
choice we can study the dynamics of the wave packet in the region where the electric 
and the magnetic fields are zero. Obviously, in this region the Lorentz force acting on 
the electron, F = e(E + v x B), is also zero. 

The dynamics of the electron motion between two infinite and parallel solenoids is 
essentially the same as the electron dynamics in the case of Aharonov-Bohm effect. In 
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Figure 4. Changes in the shape of the probability density \^\ 2 during one of the 
subsequent orbits of the wave packet. The circle represents the classical orbit. The 
animation of the dynamics can be accessed on-line |19j . 

both cases the electron moves in a field free region with no Lorentz force acting on it. 
Because of the symmetry, however, the dynamics in the case of two infinite solenoids 
consists only of the motion along one straight line between the solenoids. This avoids the 
complications of the Aharonov-Bohm dynamics where the trajectories on opposite sides 
of a single solenoid are compared. Because the properties of the dynamics are not studied 
through an interference, in the case of two solenoids the changes in the dynamics of the 
electron wave packet exclude the contribution of the nonlocal properties of the electron 
wave function dependent on the topology of the space. Another advantage of using two 
solenoids and localized wave packet is the possibility to separate the solenoids far enough 
that no parts of the wave packet penetrate the non-zero field region inside the solenoids, 
excluding this contribution to the dynamics as well. If necessary, the solenoids could be 
shielded by a potential barrier of a supercritical potential additionally preventing any 
penetration of the wave packet into non-zero field region. 

The motion of the wave packet between the solenoids, in the plane normal to 
the orientation of the solenoids, is shown in Figure [7j The wave packet was initially 
positioned away from center of the solenoids. The dynamics was studied for two initial 
momenta, p\ = 0.53 MeV/c or pi = 0.64 MeV/c. The solenoids were separated by a 
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Figure 5. The shape of the probability density l^l 2 in the vertical plane at the time 
when the wave packet crosses it. As the wave packet disperses and reaches the non- 
penetrable scalar potential Aq, the electron wave function creates a bound state in the 
direction of the magnetic field. As the wave packet rotates in the uniform magnetic 
field, the number of nodes reduces. The related animation can be accessed on-line |19j . 
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Figure 6. Vector potential of two parallel infinite solenoids with oppositely directed 
magnetic fluxes of the same magnitude. The potential is shown in the plane normal 
to the orientation of the solenoids. 



distance of 2a = 0.1 nm and the magnitude of the magnetic flux inside each of them 
was $ = 5.2 x lCr 14 Wb. 

As shown in Figure [7], and in the related animation, the wave packet moved between 
two solenoids along a straight line, dispersing in time in the direction normal to the 
direction of motion. Since, classically, there was no force acting on the electron, one 
should have expected a constant velocity of the wave packet along the straight line. 
This was not the case. The velocity of the wave packet, shown in Figure [8] as a function 
of the position, increased as the packet approached the solenoids and decreased as the 
packet left the solenoids. It is somewhat paradoxical that while there was, classically, 
no force acting on the electron, the electron acceleration, obtained strictly as a solution 
of the Dirac equation, was not zero. In the next section we will show that this paradox 
can be resolved. 



5. The dynamics of the charged particle under the conditions of the 
Aharonov-Bohm effect 

The controversy in the previous section, resulting from the solution of the Dirac 
equation, can be understood if we look again at the basic principles of classical 
mechanics. 

The motion of a particle in the absence of a force in classical dynamics is obtained 
by solving the differential equation 
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Figure 7. Six stages of the motion of the wave packet between two infinite parallel 
solenoids. Small yellow circles represent the positions and the sizes of the solenoids. 
The animation can be accessed on-line [T5] . 



where p = mv is the momentum of a particle. The solution of this equation results in a 
constant velocity of the particle. The solution of the Dirac equation shows that, even if 
the force acting on the particle is zero, the velocity of the particle changes. To resolve 
the discrepancy between classical dynamics and the Dirac equation we need to look at 
the basics of how the Dirac equation of an electron is obtained. 

In obtaining the Dirac equation from classical dynamics, the classical momentum 
p becomes a momentum operator p = —ih~V. If the electron is in a given external field 
the momentum operator is replaced by p — > —iftV + \e\A, where e is the charge of the 
electron [20]. This replacement not only is not self-evident, but must be done in the first 
order equations [20] • Going backward from the Dirac equation into classical dynamics, 
one replaces the classical momentum p — > p + \e\A. In the Lagrangian mechanics, the 
momentum p = mv + \e\A is called the generalized or canonical momentum ^T\. By 
proper use of momentum conservation [21], the motion of an electron in the absence of 
a classical force can now be obtained by solving the differential equation 

d( r + J eli) = o. (14) 
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Figure 8. The velocity, as a function of position, of a wave packet moving along a 
straight line between two infinite parallel solenoids normalized to the speed of light . 
The solenoids are positioned at x=0. 



Since A x does not explicitly depend on time, for the motion along the x-axis, Eq. (14) 
can be written as 

dp x , , dAr , , dAr dx , , dA^ 

~dt 



— \e 

In the relativistic case 

dp x _ 
dt 



dt 



— e 



dx dt 



— e 



dx 



(15) 



dv,j 



mo- 



di 



(16) 



Here m is the rest mass of the electron. Using Eqs. (14) and (16) we get the differential 
equation for the velocity of the electron 



dv x 
~dt 



\e\ 
m 



v~ 



1 V x 
1 - "I 



5 dA a 



dx 



(17) 



Due to the relativistic motion, in the system of the solenoids the velocity differential dv x 
has to be substituted with dv x /{l — v^/c 2 ). This leads to the final form of the differential 
equation 



dv x 
~dt 







—v x 


H) 


m 





dA x 
dx 



Substituting A x from the Eq. (12) for z = 0, electron motion along the mid-path 
between two solenoids, the differential equation for the velocity of the electron becomes 



dv x 
~dt 



2|e|$ 



x 



nrriQ (x 2 + a 2 ) 2 



1 



(19) 



where a is half the distance between two parallel infinite solenoids. This differential 
equation can be solved numerically. 
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Figure |9| shows a nearly perfect match between the solution of Eq. ( 19 ) and the 
solution of the Dirac equation. A small difference, a fraction of a percent, is the result of 
the numerical precision and the fact that in classical dynamics we calcute propagation 
of a particle, not a wave packet, therefore having no contribution of the localization to 
the initial condition. 

In a little digression from the main point of this paper, we can show that the 
above application of the conservation of generalized or canonical momentum is not an 
exception. One can easily show that for the translationally invariant gauges, Eqs. (|9]) 



and (10), the conservation of generalized or canonical momentum is equivalent to the 
Lorentz force. Without changing the general result, for ease of calculation, we assume 
that the electric field E = and the vector potential A is time independent. 



The conservation of generalized or canonical momentum, Eq. (14), for the dynamics 
of an electron in a uniform magnetic field defined by the translationally invariant gauges, 
can now be rewritten as 

d f t = - \e\^ = -\e\(vV)A = -\ e \(vV)(-B z i + B x k) 

= - \e\(-B v z i + B v x k) = -\e\(v x B) (20) 



Eq. (20) shows that the Lorentz force is just a consequence of generalized or canonical 
momentum conservation. While not in such an explicit form, similar results were 
obtained by Semon and Taylor [2T] . 

In a somewhat different way the same can be shown for the rotationally invariant 
gauge, Eq. (|8]). In this case we use conservation of the canonical angular momentum 

d(fx p + \e\ fx A) d(l + |e| fx A) 



dt dt 



0. (21) 



I = fx p represents classical angular momentum. Eq. (21) can be rewritten as 



dl .d(fx A) . I df -> _ dA 

— = - e : = — e — x A + r x — 

dt ' dt \ dt dt 



\e\fx(vxB)=fxF. (23) 



vx A + f x (fV)Aj (22) 

Repeating the same algebra, but using the rotationally invariant gauge of the vector 
potential A , we get 

dt 
dl ~ 

F is again the Lorentz force. With this we have shown the equivalence of the 
conservation of generalized or canonical momentum and canonical angular momentum, 
and the action of Lorentz force for all vector potential gauges of the uniform magnetic 
field. 

To conclude this section, we can look at the ongoing question if the Aharonov-Bohm 
effect is the result of a force changing the velocity of a particle passing on opposite sides 
of a single infinitely long solenoid, or there is only a quantum-mechanical phase shift 
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Figure 9. The velocity, as a function of position, of a wave packet moving along 
a straight line between two solenoids normalized to the speed of light. The circles 
represent the solution of the Dirac equation and the line represents the solution of Eq. 



(19 1 under the same conditions. The solenoids are positioned at x=0. 



[22| |23| [2U |25j [26j [27]. Quantum-mechanically such a question does not exist. The 
dynamics of the relativistic electron should be obtained only as a solution of the time- 
dependent Dirac equation. The solution of the Dirac equation shows that the velocity 
of a wave packet changes even in the region where the magnetic field is zero. Since the 
change of the velocity depends on the gradient of the vector potential, the velocity of the 
wave packet passing on the opposite sides of the solenoid will be different. Giving to the 
solution of the Dirac equation a classical picture, the requirement of the conservation 
of the canonical momentum in the region where the magnetic field is zero but the 
vector potential is not zero, also results in a velocity change of a charged particle. 
No mechanical forces need to be present in order for the velocity of the particle to 
change. From both, the solution of the Dirac equation and from the equivalence of the 
conservation of canonical momentum and the action of the Lorentz-force, one can also 
see that no work was done on the particle. The particle exits with the same energy 
regardless on which side of the solenoid it passed. 

From the studies presented in this paper, the classical picture of the Aharonov- 
Bohm effect consists of the phase shift in the wave function of the particles passing on 
opposite sides of a solenoid being attributed to time lag resulting from different evolution 
of the velocities of the particles, without any involvement of mechanical forces and any 
change in energy. A phase shift then results in an interference pattern. 
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6. Conclusion 



In conclusion, the full three-dimensional Finite Difference Time Domain (FDTD) 
method was developed to solve the Dirac equation. In this paper, the method was 
applied to the dynamics of the electron wave packet in a vector potential in the 
arrangements associated with the Aharonov-Bohm effect. The solution of the Dirac 
equation showed that the velocity of the electron wave packet changed even in the 
region where the electric and the magnetic fields were zero, and therefore no force acted 
on the electron. 

The solution of the Dirac equation agreed with the prediction of classical dynamics 
under the assumption that the dynamics were defined by the conservation of generalized 
or canonical momentum. It was shown that in the case of a uniform magnetic field, the 
conservation of generalized or canonical momentum was equivalent to the action of the 
Lorentz force. 

The studies in this paper have helped to establish a classical picture of the 
Aharonov-Bohm effect as the interference pattern resulting from the phase shift in the 
wave function of the particles passing on opposite sides of a solenoid attributed to a time 
lag resulting from different evolution of the velocities of the particles. No mechanical 
forces need to be involved and no change in energy of the particle occurs. 
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